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Abstract 

It  is  known  that  marine  sediments  can  support  both  compressional 
and  shear  waves.  However,  published  work  on  scattering  from  irregu¬ 
lar  elastic  media  has  not  exaniined  the  influence  of  shear  on  sea  bed 
scattering  in  detail.  A  perturbation  model  previously  developed  by 
the  authors  for  joint  roughness- volume  scattering  is  used  to  study  the 
effects  of  elasticity  for  three  sea  bed  types:  sedimentary  rock,  sand 
with  high  shear  speed,  and  sand  with  “normal”  shear  wave  speed. 

Both  bistatic  and  monostatic  cases  are  considered.  For  sedimentary 
rock  it  is  found  that  shear  elasticity  tends  to  increase  the  importance 
of  volume  scattering  and  decrease  the  importance  of  roughness  scat¬ 
tering  relative  to  the  fluid  case.  Shear  effects  are  shown  to  be  small 
for  sands. 

PACS  numbers:  43.30Pc,  43.60Pt,  43.30Gv,  43.30Ma 
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INTRODUCTION 

In  a  previous  paper,  the  Born  approximation  (first-order  perturbation  the¬ 
ory)  was  developed  for  sound  scattering  by  rough,  elastic  sea  beds  [13].  This 
formalism  combines  roughness  scattering  and  scattering  due  to  volume  inho¬ 
mogeneities.  An  example  was  given  showing  strong  and  complicated  elastic¬ 
ity  effects  for  roughness  and  volume  scattering  for  sea  bed  rock.  Essen  [6] 
has  used  the  Born  approximation  to  treat  sea  bed  roughness  scattering  for 
types  ranging  from  very  soft  sediment  to  basalt.  His  results  indicate  that 
shear  effects  can  be  important  for  sands  as  well  as  for  consolidated  sediments 
and  basalt.  Yang  and  Broschat  [21],  however,  have  applied  the  Born  approx¬ 
imation  and  small-slope  approximation  to  the  interface  scattering  problem 
and  give  an  example  indicating  that  shear  effects  should  not  be  important  for 
rough  sands.  Ivakin  [12]  reached  a  similar  conclusion  with  respect  to  volume 
scattering  from  sands. 

One  purpose  of  the  present  article  is  to  provide  an  evaluation  of  the 
importance  of  shear  effects  in  both  volume  and  roughness  scattering  from  sea 
bed  types  ranging  from  sand  to  rock  along  with  a  comparison  of  the  volume 
and  roughness  components  of  scattering.  Another  purpose  is  to  study  the 
dependence  of  volume  scattering  upon  the  statistics  obeyed  by  the  volume 
fluctuations,  e.g.,  correlations  between  volume  parameters. 

Use  of  the  Born  approximation  limits  the  quantitative  validity  of  the 
results  of  this  article  to  sea  beds  having  small  roughness  and  relatively  low 
levels  of  fluctuation  in  volume  parameters.  It  is  expected,  however,  that 
our  results  will  have  qualitative  validity  more  generally.  That  is,  when  the 
Born  approximation  indicates  that  shear  effects  are  large,  this  qualitative 
statement  is  likely  to  be  true  even  when  the  level  of  roughness  or  volume 
inhomogeneity  is  sufficient  to  invalidate  the  Born  approximation. 


1  SEA  BED  STATISTICAL  MODEL 

It  is  assumed  that  a  rough  interface  separates  a  lossless,  homogeneous  fluid 
(representing  sea  water)  from  an  inhomogeneous  elastic  sea  bed  medium. 
Bulk  properties  of  the  sea  bed  are  defined  by  density,  p,  and  compressional 
and  shear  wave  speeds,  Cp  and  Ct  .  Subscripts,  p  and  t,  here  and  below  denote 
longitudinal  (compressional)  and  transverse  (shear)  waves,  respectively.  In 
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the  water,  Ct  is  taken  to  be  zero. 

The  sea  bed  parameters  will  be  divided  into  two  classes:  those  represent¬ 
ing  the  average,  non-fluctuating  values  of  bulk  properties  such  as  density  and 
wave  speeds,  and  those  describing  sea  bed  random  properties.  The  former  are 
defined  by  the  seabed  type  and  will  be  assigned  values  appropriate  to  sand 
and  rock.  The  latter  are  defined  by  the  statistics  of  the  random  fluctuation 
of  bulk  parameters  and  interface  relief. 


1.1  Average  Properties 

The  average  sea  bed  parameters  of  interest  are  as  follows.  The  density  ratio. 
Op  =  p/p/,  is  the  ratio  of  the  average  sediment  mass  density  to  that  of  the 
overlying  water.  The  compressional  wave  speed  ratio,  ctp  =  Cp/c/,  is  likewise 
the  ratio  of  the  average  sediment  compressional  wave  speed,  Cp,  and  the  water 
compressional  wave  speed,  c/.  This  ratio  is  complex,  with  imaginary  part 
determining  absorption  loss.  The  shear  speed  ratio,  at  =  Ctlcj^  is  the  ratio 
of  the  average  complex  sediment  shear  speed,  Ct,  and  the  compressional  wave 
speed  in  the  water.  The  “loss  parameters”  [14,  18]  can  be  expressed  in  terms 
of  the  real  and  imaginary  parts  of  the  compressional  and  shear  speed  ratios 


Im[aQ 

Re[ac, 


(1) 


where  the  subscript  a  =  p,t  denotes  the  wave  type  in  the  seabed. 


1.2  Model  for  Roughness 

The  roughness  two-dimensional  (2-D)  spectrum  is  taken  to  be  of  the  form 


^r(K)  = 


+  K^) 


2\~tr/2 


where  K  =  K  =  yj Kl  K^.  Equation  2  assumes  that  roughness 

statistics  are  stationary  and  isotropic,  with  a  spectrum  that  obeys  a  power 
law  for  ^  Kq  .  The  parameter  Kq  is  an  inverse  correlation  scale  for 
roughness  in  the  transverse  directions.  The  roughness  spectrum  strength, 
and  the  power  exponent,  7r,  can  be  estimated,  as  shown  below,  from 
available  one-dimensional  (TD)  data. 
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The  integral  over  all  K  (including  negative  spatial  frequencies)  of  the 
spectrum  yields  the  mean-square  roughness, 

=  J  ^^;\k)cPk  (3) 

This  gives  the  following  expression  for  the  dimensionless  spectral  strength; 


B^^^  =  h^K7-\'trl2-l)ln  .  (4) 

This  requires  7r  >  2,  but  this  is  shown  below  to  be  usually  valid. 

Most  of  available  data  on  roughness  properties  are  in  the  form  of  one¬ 
dimensional  spectra,  Usually  they  can  be  well  approximated  by 

the  power  law 

j(r)L3— 

I  (5) 

The  trivial  parameter  /iq  =  1  m  is  introduced  for  convenience  to  yield  a 
dimensionless  roughness  spectrum  strength,  A^''\  Typical  measured  values 
(see,  e.g.,  [2])  are  10”®  to  10”^,  and  1  to  3. 

The  relation  between  1-D  and  2-D  spectra  is  given  by  the  integral 

$(;^(A'i)  =  J  ^i^\K)dK2. 

Then,  using  Eq.  2,  one  obtains 


2B('-) 


(A7  +  Jo  {x^  +  l)>/2  • 


Comparing  6  to  5,  one  obtains  the  relation  between  parameters  and 
7r  ^■nd  , 


Hr  —  T  1, 


(7) 


B^^^  =  hl-^^F{^r)A^’-\  (8) 

where 

""  r(i/2)r(er/2)’ 
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with  r  being  the  gamma  function.  In  particular,  one  obtains  F(l)  =  I/tt,  F(2)  = 
1/2,  F{3)  —  2/7r,  F(4)  =  3/4.  This  generalizes  the  result,  given  in  Ref.  [14] 
for  the  case  Kq  =  0.  Note  that  w  1  to  3  corresponds  to  7^  «  2  to  4. 

1.3  Model  for  Volume  Inhomogeneities 

A  description  of  the  statistics  of  volume  inhomogeneities  is  given  by  a  ma¬ 
trix  of  spatial  (3-dimensional)  cross-spectra.  These  are  Fourier  transforms 
of  cross-correlation  functions  for  the  relative  fluctuations  Cp  =  Ap/p,  Cp  = 
Acpfcp,et  =  Act/ct-  These  relative  fluctuations  will  be  collectively  desig¬ 
nated  ep  with  j3  =  p,p,t.  In  this  notation,  the  cross-spectra  are  taken  to  be 
of  the  form 


«S.(K,9)  = 


B, 


00' 


iq^o  +  q^  +  K^a'^) 


(10) 


This  equation  assumes  that  inhomogeneity  statistics  are  stationary  and  isotropic 
in  transverse  directions,  with  a  spectrum  that  obeys  a  power  law  if  -|- 
;:§>  ql.  The  parameter  qo  is  an  inverse  correlation  scale  for  inhomo¬ 
geneities  in  the  vertical  direction.  The  so-called  aspect  ratio,  a,  specifies 
the  anisotropy  of  the  inhomogeneities  and  is  the  ratio  of  horizontal  to  ver¬ 
tical  correlation  scales  [11,  24,  22,  23].  For  a  =  1,  the  inhomogeneities  are 
isotropic  on  the  average  and  have  a  spherical  form  (referring  to  surfaces  of 
constant  correlation  in  spatial  lag  space).  For  a  »  1,  the  inhomogeneities  are 
strongly  anisotropic  and  constant  correlation  surfaces  are  thin  oblate  ellip¬ 
soids.  Generally,  in  the  case  of  transverse  anisotropy,  one  more  aspect  ratio 
must  be  introduced  corresponding  to  the  third  correlation  scale  [4].  Angular 
coordinates  for  orientation  of  the  inhomogeneity  ellipsoids  in  space  can  also 
be  arbitrarily  introduced  [22,  23]. 

The  three-dimensional  integrals  analogous  to  Eq.  3  are  the  elements  of 
the  covariance  matrix  for  volume  fluctuations 


<  £/)£S.  >=  /  4'S.(K,  q)d^K  dq.  (11) 

These  integrals  are  infinite  if  7„  <  3  but  can  be  made  finite  for  all  jy  if  a 
high  frequency  cutoff  is  introduced.  Accordingly,  it  will  be  assumed  that 

=  0  (12) 
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for 


K'^a^  +  q^  >  {qolef,  e  <  1.  (13) 

Introduction  of  the  small  dimensionless  cutoff  parameter,  e,  provides 
smooth  behavior  of  the  corresponding  correlation  functions  for  small  spa¬ 
tial  lags  and  also  yields  finite  variance  for  the  inhomogeneity  fluctuations  for 
any  7„.  Obviously,  this  is  important  as  experimental  data  have  finite  vari¬ 
ance  and  have  finite  spatial  resolution  with  respect  to  both  large  and  small 
scales.  The  “small”  spatial  lags  for  which  the  cutoff  has  effect  have  scales 
comparable  to,  or  smaller  than  the  components  of  the  vector  (o,  a,  l)e/go- 
The  matrix  elements  Bpf}/  can  be  determined  from  Equation  11  (see  Ap¬ 
pendix).  Otherwise,  available  geoacoustic  data  can  be  used  for  the  estimation 
of  the  spectrum  parameters.  In  the  case  of  strongly  correlated  or  anticor¬ 
related  fluctuations  one  can  consider  any  one  of  the  three  fluctuations,  say 
density,  as  defining  the  other  two.  This  parameter  will  be  called  the  “refer¬ 
ence”  parameter.  In  this  case,  put 


ep  =  rfipCp,  (14) 

where  the  corresponding  speed  fluctuation  ratios  (with  respect  to  fluctuations 
of  the  density)  are 

where  rpp  can  be  taken  from  compilations  of  geoacoustic  data,  for  example 
Hamilton’s  [7,  9,  10,  8]  and  others  [17,  20].  In  this  case 

Bpp'  —  t'pp  t pi p  Bpp,  (fb) 

for  all  combinations  of  (3  and  j3\  with  r^p  =  1. 

The  parameters  Bpp  and  ')y  can  be  also  obtained  from  available  data. 
Most  data  on  fluctuations  of  bulk  properties  provide  only  one-dimensional 
(vertical)  spectra,  $i’^^(9)  .  Usually  they  can  be  well  approximated  by  the 
power  law 


=  — f — ,  9»go.  (17) 

Again,  the  parameter  /iq  =  1  ni  is  used  to  obtain  a  dimensionless  spectrum 
strength,  Typical  measured  values  (see,  e.g.,  [1,  3, 16,  22,  24])  are  « 
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10  to  10  ,  fa  0.5  to  2.  Typical  estimates  obtained  from  inversions  of 

geoacoustical  data  [11,  22,  24]  are  a  «  1  to  10. 

The  relation  between  1-D  and  3-D  spectra,  generalized  in  Ref.  [24]  for 
anisotropic  inhomogeneities  with  aspect  ratio  a,  is  given  by  the  expression 


where 


27Tq'  dq'  ’ 


=  sjq^  +  K^a\ 

Then,  using  Eq.  17,  one  obtains  the  relation  between  parameters  and 
'jy  and 


Iv  —  A  2, 

& 

=  hl~^'’A^a'^(yl{2n). 

2  NUMERICAL  EXAMPLES 

To  illustrate  the  effects  of  shear  elasticity,  we  consider  scattering  from  sea 
beds  of  varying  composition,  ranging  from  sand  to  solid  rock.  Finer-grained 
silts  and  clays  are  not  considered,  as  their  low  shear  speeds  [6,  19]  result  in 
negligible  shear  effects.  The  parameters  Cp,  ct,  and  p  that  characterize  these 
examples  were  chosen  primarily  to  illustrate  certain  points  rather  than  to 
provide  extensive  coverage  of  a  range  of  realistic  sea  bed  types.  Some  of  the 
chosen  examples  are  fairly  realistic,  others  are  not.  The  parameters  spec¬ 
ified  below  were  used  in  computing  angular  dependencies  of  the  scattering 
strength.  With  regard  to  roughness  scattering,  some  of  the  effects  seen  in 
the  present  example  are  visible  in  other  published  work  [5,  21,  6].  The  intent 
here  is  to  compare  and  contrast  elasticity  effects  as  they  appear  in  roughness 
and  volume  scattering. 
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2.1  Scattering  Strength 

In  this  section,  the  expressions  given  in  [13]  are  used  for  calculation  of  scat¬ 
tering  strength,  10  log  c  ,  where  a  is  the  scattering  cross  section  per  unit 
area  of  the  sea  bed  and  characterizes  the  frequency-  angular  distribution  of 
the  mean  intensity  of  field  fluctuations  in  the  Fraunhofer  zone  (or  far  zone, 
relative  to  a  scattering  surface  or  volume).  The  scattering  cross  section  is 
assumed  to  be  the  sum  of  contributions  from  interface  roughness  and  volume 
inhomogeneity  in  the  form 

where  <7^  and  ay  are  components  due  to  roughness  and  volume  scattering, 
respectively,  and  K,  are  horizontal  (transverse)  components  of  the  scat¬ 
tered  and  incident  wave  vectors,  correspondingly.  Vectors  and  K,-  are 
expressed  through  the  essential  angular  coordinates  of  the  problem 

=  {kf  cos  0s,i  cos  kf  cos  9s,i  sin  4>s,i)  •  (18) 

These  are:  grazing  angles,  Oi  and  ds,  for  the  incident  and  scattered  acoustic 
waves,  respectively,  and  the  azimuths,  (t>s  and  (^,  for  the  incident  and  scattered 
waves  (see  Figure  1).  Frequency  enters  through  the  acoustic  wave  number 
in  the  water  kj  =  ojjcj.  Roughness  and  volume  scattering  cross  sections, 
cTr  and  <7„,  are  expressed  in  terms  of  the  average  parameters  of  the  seabed, 
ap,ap,at,  and  its  statistical  parameters,  the  two-dimensional  spatial  spectra 
of  roughness,  $*'')(K),  and  the  three-dimensional  cross-spectra  of  different 
volume  fluctuations,  $^'jg,(K,5),  where  K  =  —  K,  and  q  is  the  vertical 

component  of  the  spectral  argument  [13]. 

2.2  Input  parameters 

Table  1  lists  the  average  properties  of  the  three  sea  bed  types  chosen  for 
illustration  of  the  effects  of  elasticity  on  acoustic  scattering.  These  sea  bed 
types  include  a  strongly  elastic  case  (sedimentary  rock),  a  moderately  elastic 
case  (sand  with  unusually  high  shear  speed,  “shear”  sand),  and  a  weakly 
elastic  case  (sand  with  a  more  typical,  “moderate”,  shear  speed). 

Two  different  sets  of  random  parameters  are  used  for  sands  and  rock. 
The  roughness  and  volume  parameters  are  specified  below. 
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The  roughness  spectral  parameters  were  assigned  the  values  used  in  Ref 
[13]:  7r  =  4,  Ao  =  and  5^’')  =  (2/7r)  x  with  =  10“®,  with 

the  factor  {2/n)  resulting  from  the  relation  between  2-D  and  1-D  roughness 
spectra. 

The  reference  (density)  inhomogeneity  spectrum  was  taken  to  be  of  the 
form  given  in  Eq.  10  with  7v  =  3,  qo  =  10“^m“\  =  (27r)“^a^  x  and 

=  10“^.  Again,  these  parameters  were  used  in  Ref.  [13],  except,  here, 
the  aspect  ratio  is  not  taken  to  be  unity.  The  parameter  contains  a 
factor  a? I{2'k)  resulting  from  the  relation  between  3-D  and  1-D  spectra.  The 
corresponding  one-dimensional  (vertical)  spectrum  is  simply  A^'^'^jK,  which 
gives  a  reasonable  fit  to  the  available  data  [24]. 

In  the  “correlated”  case,  the  cross  spectrum  matrix  elements  were  cal¬ 
culated  using  Eq.  16  with  the  ratios  r^p  estimated  from  published  relations 
between  geoacoustic  parameters  [10,  17]  and  listed  in  Table  2  for  the  sea  bed 
types  chosen  for  illustration. 

In  another  case  to  be  considered,  the  elements  of  the  correlation  matrix 
are  taken  to  be  of  the  form 

=  (19) 

Because  of  the  Kroneker  delta  used  in  Eq.  19,  all  the  nondiagonal  elements 
of  the  correlation  matrix  vanish.  In  this  “uncorrelated”  case,  the  volume 
component  of  the  scattering  cross  section,  according  to  [13],  will  be  simply 
the  sum  of  three  components 

=  (20) 

where  correspond  to  contributions  of  fluctuations  of  the  three  different 
bulk  parameters  {(3  =  p,Cp,Ct). 

Calculations  are  presented  both  for  the  sum  and  for  each  of  the  terms  in 
Eq.  20.  The  former  permits  an  evaluation  of  the  role  of  cross  correlations 
between  different  volume  fluctuations,  or  nondiagonal  elements  of  the  corre¬ 
lation  matrix,  by  comparison  of  the  two  extremes:  the  strongly  correlated 
case  and  the  uncorrelated  case.  The  latter  permits  a  comparison  of  the  the 
role  of  different  fluctuations. 
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3  DISCUSSION 

The  first  example  to  be  considered  uses  parameters  appropriate  to  sedimen¬ 
tary  rock  [6].  As  the  set  of  figures  used  in  this  and  following  examples  is  rather 
involved,  an  explanation  of  the  arrangement  of  the  figures  will  be  given  before 
discussing  the  results.  Figure  2  shows  the  monostatic  or  backscattering  case, 
giving  the  dependence  of  scattering  strength  on  backscattering  grazing  angle. 
Figure  3  shows  the  dependence  of  bistatic  scattering  strength  upon  scattered 
grazing  angle  with  the  incident  grazing  angle  fixed  at  45°  and  the  difference 
between  scattered  and  incident  azimuth  fixed  at  180°.  As  the  scattered  graz¬ 
ing  angle  varies  from  0°  to  180°,  it  passes  through  the  backscatter  direction 
(45°)  and  the  forward,  or  “specular”,  direction  (135°).  Figure  4  shows  the 
dependence  of  bistatic  scattering  strength  upon  the  difference  between  scat¬ 
tered  and  incident  azimuthal  angles  at  fixed  incidence  and  scattered  grazing 
angles,  both  taken  to  be  45°.  All  three  figures  compare  the  elastic  and  fluid 
cases.  In  the  fluid  case,  the  shear  wave  speed  is  set  to  zero  and  all  other  pa¬ 
rameters  are  kept  the  same  as  in  the  elastic  case.  Each  of  the  figures  is  divided 
into  parts  (a)  and  (b).  Part  (a)  compares  three  separate  cases:  roughness 
scattering,  correlated  volume  scattering  and  uncorrelated  volume  scattering. 
Part  (b)  compares  separate  volume  scattering  mechanisms,  fluctuations  in 
density,  compressional  wave  speed,  and  shear  wave  speed. 

For  sedimentary  rock,  the  effects  of  shear  elasticity  on  both  roughness  and 
volume  scattering  are  very  strong,  as  illustrated  by  Figs.  2-4.  Noteworthy 
effects  are: 

-  Figures  2a  and  3a  show  that  roughness  scattering  strength  exhibits 
dips  at  grazing  angles  immediately  below  the  critical  angles  for  shear  and 
compressional  waves  (39.7°  and  64.2°,  respectively). 

-  Shear  elasticity  decreases  roughness  scattering  strength  compared  to  the 
fluid  case  for  grazing  angles  less  than  the  shear  critical  angle  and  near  the 
forward  direction  (Figs.  2a,  3a,  and  4a). 

-  Elasticity  greatly  enhances  scattering  strength  for  volume  scattering 
in  the  directions  where  it  decreases  roughness  scattering  (Figs.  2-4).  The 
effects  for  volume  scattering  are  much  larger  than  the  effects  for  roughness 
scattering. 

-  For  grazing  angles  below  the  compressional  critical  angle,  scattering  due 
to  density  fluctuations  has  nearly  the  same  angular  dependence  as  scattering 
due  to  shear  speed  fluctuations  (Figs.  2b  and  3b). 
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-  The  effect  of  cross-correlation  for  different  perturbations  can  be  large, 
as  it  is  noticeable  in  Fig.  3a  from  the  difference  between  correlated  and 
uncorrelated  cases  in  the  range  60°  -  120°. 

The  influence  of  shear  elasticity  in  surficial  sands  is  expected  to  be  rather 
small  as  they  typically  have  shear  speeds  much  lower  than  the  water  com- 
pressional  wave  speed  (see,  e.g.,  [10,  20]).  For  illustrative  purposes,  therefore, 
an  example  will  first  be  considered  in  which  shear  effects  are  exaggerated. 
This  example  is  used  by  Essen  [6]  and  is  distinguished  by  a  shear  speed  of 
480  m/s.  This  case  of  sand  with  high  shear  speed,  the  “shear”  sand  case,  is 
illustrated  in  Figs.  5-7.  These  figures  show  the  following  effects: 

-  Both  roughness  and  volume  backscattering  (Fig.  5a)  are  reduced  relative 
to  the  fluid  case  for  grazing  angles  smaller  than  the  compressional  critical 
angle  (33.6°).  Compared  to  the  sedimentary  rock  example,  this  is  the  same 
behavior  regarding  roughness  but  the  opposite  regarding  volume  scattering. 

-  The  relative  contributions  of  volume  and  roughness  scattering  (Figs. 
5a,  6a,  7a)  are  very  dependent  on  scattering  angle.  Roughness  scattering  is 
dominant  near  the  specular  direction,  but  volume  scattering  can  be  important 
for  incident  or  scattering  angles  greater  than  the  compressional  critical  angle. 

-  The  effect  of  cross-correlation  for  different  perturbations  can  be  large,  as 
it  is  noticable  from  the  difference  between  correlated  and  uncorrelated  cases 
in  Fig.  6a  (in  the  range  80°  to  180°)  and  in  Fig.  7a  (-100°  to  100°).  Note 
that  this  difference  is  negligible  in  the  fluid  case. 

-  Volume  scattering  in  backward  directions  for  sands  (Figs.  5b,  6b,  7b) 
is  mostly  due  to  density  fluctuations,  while  both  compressional  and  shear 
wave  speed  fluctuations  are  unimportant.  But,  at  the  same  time,  even  weak 
fluctuations,  of  compressional  speed  (an  order  of  magnitude  smaller  than  den¬ 
sity  fluctuations)  can  be  important  near  the  forward  direction  (Fig.  6b). 
Shear  fluctuations  also  provide  significant  scattering  near  the  forward  direc¬ 
tion  (Figs.  6b,  7b). 

The  last  example,  with  scattering  strength  shown  in  Figs.  8  and  9,  is 
a  sand  with  shear  speed  of  about  200  m/s,  similar  to  typical  measured  val¬ 
ues  [10,  19].  Here  again,  the  relative  importance  of  roughness  and  volume 
scattering  is  dependent  upon  angle.  Density  fluctuations  are  generally  more 
important  than  compressional  speed  fluctuations  in  the  backscattering  case 
(Fig.  8)  with  relatively  reduced  contribution  in  forward  directions  (Fig.  9). 
These  effects  have  little  to  do  with  shear  elasticity,  however.  In  fact,  the 
effects  of  shear  elasticity  are  negligible  for  both  roughness  and  volume  scat- 
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tering.  Scattering  by  shear  speed  fluctuations  in  forward  and  near  specular 
directions  is  about  10  dB  lower  than  scattering  by  density  and  compressional 
speed  fluctuations  (Fig.  9)  and  absolutely  negligible  in  the  backscattering 
case  (the  shear  component  is  not  seen  in  Fig.  9  as  its  level  is  too  low).  Thus, 
even  though  sand  is  definitely  a  shear-supporting  medium  in  a  geotechnical 
sense,  it  can  be  treated  by  a  simple  fluid  acoustic  model. 

The  effects  of  elasticity  on  roughness  and  volume  scattering  for  rock  can 
be  summarized  by  saying  that  elasticity  “softens”  the  interface,  decreasing 
the  acoustic  contrast  with  water  with  a  corresponding  decrease  in  reflection 
and  roughness  scattering  and  a  corresponding  increase  in  transmission  and 
volume  scattering.  This  view  is  substantiated  by  Fig.  10,  which  compares 
the  magnitudes  of  reflection  and  transmission  coefficients  for  the  fluid  and 
elastic  cases.  In  the  fluid  case,  the  reflection  coefficient  is  very  nearly  equal  to 
unity  for  grazing  angles  below  the  compressional  critical  angle.  The  inclusion 
of  lossy  shear  greatly  alters  this  picture,  with  the  energy  reflection  coefficient 
dropping  to  very  small  values  for  grazing  angles  in  the  vicinity  of  35°.  This 
causes  the  decrease  noted  in  roughness  scattering.  At  the  same  time.  Fig. 
10  shows  that  the  transmission  coefficient  for  pressure  waves  is  greatly  in¬ 
creased  when  elastic  effects  are  included.  The  conversion  to  shear  energy 
is  also  efficient,  and  both  effects  enhance  volume  scattering.  The  opposite 
effect  for  scattering  in  high-shear-speed  sand,  mentioned  above,  is  explained 
analogously  by  Fig.  11. 

For  the  sedimentary  rock  example,  when  the  incident  and  scattered  graz¬ 
ing  angles  are  less  than  the  compressional  critical  angle,  only  incident  and 
scattered  vertically  polarized  shear  waves  are  important.  Thus,  there  is  no 
interference  between  different  scattering  channels;  only  shear-to-shear  scat¬ 
tering  is  important.  This  causes  the  similarity  observed  in  the  angular  depen¬ 
dences  of  scattering  by  density  and  shear  speed  fluctuations.  The  scattering 
levels  due  to  these  two  mechanisms  are  different  only  because  different  fluc¬ 
tuation  levels  are  assumed  (Table  2). 

4  CONCLUSIONS 

First-order  perturbation  theory  results  obtained  in  [13]  have  been  used  to 
study  the  effects  of  shear  elasticity  on  acoustic  scattering  by  sea  beds  of  dif¬ 
ferent  types.  While  the  calculations  presented  here  are  quantitatively  correct 
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only  in  a  region  of  small  roughness  and  small  volume  fluctuations,  it  is  likely 
that  the  general  conclusions  drawn  from  these  examples  will  be  qualitatively 
valid  outside  this  region. 

With  respect  to  scattering  by  rock,  shear  elasticity  weakens  roughness 
scattering  and  strengthens  volume  scattering.  Surprisingly,  this  suggests 
that  volume  scattering  may  sometimes  dominate,  or  at  least  compete  with, 
roughness  scattering,  in  sharp  contrast  to  the  predictions  of  the  fluid  model. 
The  other  general  conclusion  resulting  from  these  examples,  is  that  sand  can 
be  treated  as  an  acoustic  fluid.  This  reinforces  the  results  presented  by  Yang 
and  Broschat  [21]. 
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Appendix:  Volume  Cross  Spectrum  Matrix 

The  matrix  elements,  can  be  determined  from  Equation  11  which  gives 

where 


<^/3  =  V<  VH?  > 

are  the  standard  deviations  of  the  three  different  fluctuation  types,  and 


R/3/3'  = 


<  > 

(7  per  pi 


are  their  cross-correlation  coefficients.  In  the  case  of  strongly  correlated  or 
anticorrelated  fluctuations  one  can  put  Rppi  =  ±1  respectively,  and  in  the 
“uncorrelated”  case,  Rppi  =  Sppi.  The  parameter  is  a  dimensionless 

normalizing  factor 


/•l/e 

Q(l,e)  =  2 

Jo 


(1  -h  x2)T'/2- 


dx. 


For  7„  =  3,  one  obtains 
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while,  for  7^,  <  3, 


(5  «  2  ln(l/£)  , 


Q«2e-P-”>/(3-7.). 

and,  for  7^  >  3, 

Q«^(7v,0)  =  r(l/2)r(7./2-3/2)/r(7./2-l)  ,  (21) 

For  7v  =  4,  in  particular,  Q  =  n.  Expression  21  is  exact  if  e  =  0. 
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Table  1:  Parameters  defining  average,  non-fluctuating  bulk  properties  of  the 
illustrative  sea  bed  types. _ _ _ 


Descriptive 

Name 

Density 

Ratio 

Compressional 

Speed 

Ratio 

Shear 

Speed 

Ratio 

Reference 

m 

at 

“Normal”  Sand 

2.0 

1.20  -  0.02i 

0.133  -  O.Oli 

Hamilton  (1980) 

“Shear”  Sand 

2.0 

1.20  -  0.005i 

0.32  -  0.07i 

Essen  (1994) 

Sedimentary  Rock 

2.5 

2.30  -  0.004i 

1.30  -  O.lli 

Essen  (1994) 

Table  2:  Parameters  defining  random  properties  of  the  illustrative  sea  bed 
types. _ 


Volume 

Volume 

Volume 

Compress. 

Shear 

Rough. 

Rough. 

Sea  Bed 

Spectr. 

Spectr. 

Aspect 

Speed  Fluct. 

Speed  Fluct. 

Spectr. 

Spectr. 

Type 

Strength 

Expon. 

Ratio 

Ratio 

Ratio 

Strength 

Expon. 

iv 

a 

^tp 

^(r) 

Sand 

10-^ 

1 

3 

0.1 

2.0 

10"® 

3 

Rock 

10-^ 

1 

3 

3.0 

4.0 

10-^ 

3 
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FIGURES 


Figure  1.  Geometry  relevant  to  bistatic  scattering  from  a  random  sea  bed. 

Figure  2.  Backscattering  strength  as  a  function  of  grazing  angle  for  the 
sedimentary  rock  example.  The  fluid  case  uses  identical  parameters,  except 
the  shear  wave  speed  is  set  to  zero.  The  abscissa  is  6s  =  6i  with  4>s—4>i  =  180°. 
With  the  parameters  used,  there  is  no  frequency  dependence.  Here  and 
below:  (a)  roughness  scattering,  total  volume  perfectly  correlated  and  total 
volume  uncorrelated  scattering,  (b)  scattering  from  different  types  of  volume 
inhomogeneities. 

Figure  3.  Bistatic  scattering  strength  as  a  function  of  scattered  grazing 
angle  for  the  sedimentary  rock  example.  The  grazing  angle  of  the  incident 
wave  is  45°.  The  abscissa  is  Bg  with  —  (f>i  —  180°.  All  parameters  are  the 
same  as  for  Fig.  2. 

Figure  4.  Bistatic  scattering  strength  as  a  function  of  scattered  azimuth, 
0s  —  0,-,  for  the  sedimentary  rock  example.  The  grazing  angles  of  the  incident 
and  scattered  waves  are  Bs  =  6i  =  45°.  All  parameters  are  the  same  as  for 
Fig.  2. 

Figure  5.  Backcattering  strength  as  a  function  of  grazing  angle  for  the 
high-shear-speed  sand  example.  The  fluid  case  uses  identical  parameters, 
except  the  shear  wave  speed  is  set  to  zero. .  The  abscissa  is  Bs  =  Bi  with 
4>s  —  4>i  =  180°.  With  the  parameters  used,  there  is  no  frequency  dependence. 

Figure  6.  Bistatic  scattering  strength  as  a  function  of  scattered  grazing 
angle  for  the  high-shear-speed  sand  example.  The  abscissa  is  Bs  with  4>s—<pi  = 
180°.  All  parameters  are  the  same  as  for  Fig.  5. 

Figure  7.  Bistatic  scattering  strength  as  a  function  of  scattered  azimuth, 
0s  —  0i,  for  the  high-shear-speed  sand  example.  The  grazing  angles  of  the 
incident  and  scattered  waves  are  Bs  =  Bi  =  45°.  All  parameters  are  the  same 
as  for  Fig.  5. 

Figure  8-  Backcattering  strength  as  a  function  of  grazing  angle  for  the 
“normal”  sand  example.  The  fluid  case  uses  identical  parameters,  except  the 
shear  wave  speed  is  set  to  zero.  The  abscissa  is  Bg  =  Bi  with  0*  —  0,  =  180°. 
With  the  parameters  used,  there  is  no  frequency  dependence. 

Figure  9.  Bistatic  scattering  strength  as  a  function  of  scattered  grazing 
angle  for  the  “normal”  sand  example.  Incident  grazing  angle  Bi  =  45°.  The 
abscissa  is  Bg  with  0^  —  0j  =  180°.  All  parameters  are  the  same  as  for  Fig.  8. 

Figure  10.  Transformation  (reflection  and  transmission)  coefficients  for 
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the  sedimentary  rock  example. 

Figure  11.  Transformation  (reflection  and  transmission)  coefficients  for 
the  high-speed  sand  example. 
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